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Angular correlations in dense solutions and melts of flexible polymer chains are investigated 
with respect to the distance r between the bonds by comparing quantitative predictions of 
perturbation calculations with numerical data obtained by Monte Carlo simulation of the 
bond-fluctuation model. We consider both monodisperse systems and grand-canonical (Flory- 
distributed) equilibrium polymers. Density effects are discussed as well as flnite chain length 
corrections. The intrachain bond-bond correlation function P{r) is shown to decay as P{r) ~ 
l/r"^ for ^ <C r <C r* with ^ being the screening length of the density fluctuations and r* ^ N^^^ 
a novel length scale increasing slowly with (mean) chain length N. 



1 Introduction 

Background. It is generally assumed that large scale correlations are screened in dense so- 
lutions of flexible polymers beyond the local correlation length ^ characterizing the decay of 
the density fluctuations.^"^ One consequence of this screening hypothesis is that orientational 
correlations between two bonds Zq and li on the same chain should vanish rapidly for distances 
r ^ ^ and for corresponding curvilinear distances s = ni — uq ^ g with g being the number of 
monomers spanning the correlation length. See Figure [T] for a sketch of the notations used in 
this paper with n denoting the monomer index, the number of monomers per chain, R{N) 
the root-mean-square end-to-end distance, p the monomer number density, I the root-mean- 
square bond length, b = liniN^oo RiN)/N^/'^ the effective bond length and Cqo = 
dimensionless chain stiffness parameter.^' ^ 

Surprisingly, recent numerical studies^"^ have demonstrated the power-law decay of the in- 
trachain bond-bond correlation function P{s) = {Iq ■ li) /P, averaged over bond pairs of same 
curvilinear distance s, as a function of s 

P{s) ^ cps-'^ for <C s < (1) 

with an exponent w = 3/2 rather than the exponential cut-off expected from Flory's ideality 
hypothesis.^'' (The amplitude cp is given in eq [3] below.) This result has been rationalized 
by means of scaling arguments and perturbation calculations which demonstrate the systematic 
swelling of the chain segments.^' ''' ^^^^^ g[gi gf ^j-^g calculation is that the effective interactions 
between the monomers of a chain are only partially screened and represented (to leading order) 
by an effective potential in momentum space 

Hk) ~ ^ for e « l/k « R{N) (2) 
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(1) power law (2) finite-N correction 



Figure 1: Sketch of the considered problem in real space (a) and of the first-order perturbation 
interaction diagrams in reciprocal space (b). (a) The angular correlations are characterized by 
the bond-bond correlation function P{r) = {Iq ■ li) /t^ averaged over all pairs of bonds Zq ^'iid 
li of a chain of same distance r = \r\ and normalized by the mean-square bond length (b) 
The bold vertical arrows represent the bond vectors (eg I1U|) . dashed lines the effective monomer 
interactions ^(A:) and thin lines the Fourier-Laplace transformed Gaussian propagators G{k,t) 
with the Laplace variable t being conjugated to the curvilinear distance s. The inserted wavevec- 
tor q is conjugated to the distance r between both bonds in real space. Angular correlations 
in ring polymers are described by diagram (0). The asymptotic behavior of long linear chains 
(eq[5]) is obtained from diagram (1). The last diagram (2) describes the finite-size corrections 
relevant for large distances r r*{N) N^^^. 



increasing quadratically with wavevector /c.^'^^ The detailed calculation yields the power-law 
amplitude 

3 c, 

which is very close to the empirical values for all simulation models tested.^ 



CP = \/^— (3) 



Aim and key results of this study. Since the power-law decay of P{s) resembles the return 
probability of a random walk in three dimensions it is tempting^^ to attribute the observed 
effect to "self-kicks" of the chain involving the correlated bonds themselves (or their immediate 
neighbors). Accordingly, the bond-bond correlation function should reveal a 5(r)-correlation if 
sampled as a function of the distance r = |r| between bond pairs. This interpretation turns out 
to be incorrect, however, and we will show that the power law in s simply translates as 

s-^ ^ (rjhyl^ (4) 

with a Flory exponent = 1/2 for (to leading order) Gaussian chains. More specifically, it will 
be demonstrated by means of analytical theory and Monte Carlo simulation that 

P(r) ^ Pa(r) = ^ for g < r < r*\N) (5) 

as suggested by eq HI i.e. the angular correlations are genuinely long-ranged. The index "a" 
emphasizes that P^^r^ is the predicted asymptotic behavior for long chains on scales where the 
system behaves as an incompressible solution. For systems with finite compressibility (^ ^ b) 
the power law generalizes naturally to 

Pir) = P.iO f{n) (6) 
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in terms of the reduced distance u = r/^ and a scaling function f{u) =^ for 1 ^ u. 

Interestingly, the upper cut-off r*{N) of eqOis found to increase rather slowly with chain length 

r*(iV) 6iV^/^ < i?(iV) w (7) 

The simulation of computationally challenging chain lengths thus is required to demonstrate 
numerically the predicted power-law decay of P{r). 

Outline. The one-loop perturbation calculation leading to the above results follows again the 
seminal work by Edwards.^ See Figure (U^b) for a sketch of the computed interaction graphs. 
This calculation will be discussed first (Section [2|). Computational methods and parameters are 
summarized in section [3l We present then in section H] our numerical results obtained for systems 
containing either monodisperse polymers or Flory size-distributed equilibrium polymers. ^^'^^ By 
varying the density we scan the screening length ^ over two orders of magnitude.^ This puts us 
into a position to test the general scaling relation, eqEl for systems with finite compressibility. 
Focusing first on the properties of asymptotically long chains we consider finally the strong 
finite chain-size effects predicted by eq [71 A synopsis of our results is given in section [5] where 
we suggest possible avenues for future studies. 

2 Perturbation calculation 

General remarks. We remind^ that the first-order perturbation calculation of a quantity A 
under a perturbation U (to be specified below) generally reads (A) ~ (•^)o + (^)o (-^)o ~ (^-^)o- 
Averages performed over an unperturbed reference system of Gaussian chains of effective bond 
length b are denoted (. . .)q. In this study we have to average the observable A = Iq ■ li/t^ over 
all intrachain bond pairs at a given distance r = |r|. For linear chains {A)q = by construction. 
Thus we only have to compute the average 

{A) « - {UA), (8) 

which simplifies considerably the task compared to the perturbation calculation of the mean- 
square segment size R'^{s) presented in refs 5,7,13. We remind that for closed cycles the ring 
closure implies long range angular correlations even for Gaussian chains, hence, 

{A) « (^)o / (9) 

to leading order. 

We suppose first that the chains are infinite and on local scale perfectly fiexible (coo = {b/l)'^ = 
1).^'' We begin by formulating the problem in reciprocal space and demonstrate then the power- 
law asymptote, eq[5j Finite chain-size effects are first discussed for Flory distributed chains and 
then by means of inverse Laplace transformation for monodisperse melts. Finally, it is shown 
how our results can be reformulated for semifiexible chains (coo > 1). 

Reciprocal space description of flexible chains. The calculation of eq [8] and [9] is most 
readily performed in reciprocal space as sketched in Figure mb). The Fourier transform of a 
function f{r) is denoted f{k) = J^[f{r)] = f dr/(r)e~''°"^, the Laplace transform of a function 
f{s) is written f{t) = C[f{s)] = /(s)e~** with t being the Laplace variable conjugated to 
the arc- length s. 

Bond vectors and observable. The vertical bold arrows in Figure represent the Fourier 
transformed bond vectors 



3 



with B{1) and B{k) = J^[B{1)] being the Gaussian distribution function of the bond vector 
in real and reciprocal space, respectively. We have used here that B{k) can be expanded at 
low momentum. The wavevectors conjugated to the bonds Iq and Zi are denoted qo and qi, 
respectively. The Fourier transform of the observable A reads, hence, 

f 

^(9o,qi) = -^Qo • gi- (11) 

Effective potential v(k). The perturbation potential U{r) in real space is supposed to be the 
pairwise sum of the effective monomer interactions v{r) of all pairs of monomers of the same 
chain. To obtain v{r) one labels a few chains. The interaction between labeled monomers is 
(partially) screened by the background of unlabeled monomers. It has been shown^'^^ that 
within linear response this corresponds to an effective potential v{k) = ^[v{r)] with 

^ ^ +F{k). (12) 



v{k)p vp 

This potential is represented by the dashed lines in the diagrams. The bare excluded volume v 
indicated in the first term of eq [12] characterizes the short-range repulsion between monomers. 
Thermodynamic consistency requires^' that v is proportional to the inverse of the measured 
compressibility of the solution 

v = -^^ {a/e (13) 

where we have introduced a convenient monomeric length a = hj and defined the screening 
length following, e.g., eq 5.38 of ref 2. Please note that g can be determined experimentally 
or in a computer simulation from the low-wavevector limit of the total monomer structure 
factor and, due to this operational definition, g is called "dimensionless compressibility".^ The 
single chain form factor F{]i) represents the interaction between two monomers caused by the 
chain connectivity.^ For Gaussian chains of finite length the form factor is given by Debye's 
function.^ For infinite chains we have, hence, i^(fc) = 2/(a/c)^. For scales larger than the 
screening length (^fc <C 1) the finite compressibility (indicated by the first term in eq ll2p becomes 
negligible and the solution behaves for all densities as an incompressible melt. Equation [T2] 
reduces thus to the scale- free interaction potential already mentioned (eq[2|). 

Fourier- Laplace transform of the propagator. We remind that the Fourier transform of the 
Gaussian propagator^ may be written G{k,s) = exp(— s(a/c)^) with s denoting the curvilinear 
distance between two monomers of the chain. An arrow along a chain contour (thin lines) 
corresponds to the Fourier-Laplace transformed Gaussian propagator G{k,t) = C[G{k,s)] = 
l/{{ak)'^ + t) of wavevector k as specified in the diagrams. We need to average over all bond 
pairs at a given distance irrespective of their curvilinear distance s and we have thus to sum below 
over all possible s. For infinite chains this corresponds to setting t = for the corresponding 
Laplace variable and we shall often use the summed up Gaussian propagator G{k) = G{k, t = 
0) = (a/c)"^, i.e. the Fourier transform of the density G{r) = l/(47rra^) around a reference 
monomer of all monomers belonging to the same infinite Gaussian chain. ^ Hence, 

for infinite chains on large scales {^k <C 1). 

Interaction diagrams. The momentum q inserted in the interaction diagrams is conjugated 
to the distance r between both bonds. Momentum is a conserved quantity flowing from one 
correlated point to the other. If the momentum flows in the opposite direction of a bond (as it 
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is the case for the second bond Zi) the wavevector comes with a negative sign in eq[TTJ The first 
two diagrams are thus given by the convolution integrals 



Io{q) = [ G{qo)A{qo,-qi)G{qi) (15) 

h{q) = [ G{qo)A{qo,-qi)Giqi) 

Jqo+qi=q 

x[-v{qi)]G{qi) (16) 

The integrals Io{q) and Ii{q) describe the bond-bond correlation functions of closed cycles and 
linear chains, respectively. The perturbation calculation of linear chains using eq [8] implies a 
minus sign. This sign is indicated in front of the effective interaction v in the last equation. 
Using eq [11] and assuming eq O the integrals are considerably simplified 

hiq) = -^Io{q) (17) 



2p 
18p 

and the inverse Fourier transforms are thus 



qoG{qo) ■ qiG{qi) (18) 



qo+qi=q 



h{r) = -^Io{r) (19) 
= ^pidrG{r)f. (20) 

Sum rule for closed cycles and linear chains. Up to a constant prefactor the integrals Iq 
and Ji are thus equal on large scales {r ^ This remarkable result needs further discussion 
beyond the technical notions set up in the preceding paragraph. 

Closed cycles. The bond-bond correlation function of closed rings Po{r) is directly obtained 
from Io{r) after normalization with the density G{r) 

P„(o = #i = -fiV. (21) 



G2(r) \3r 

The reason for the normalization factor G'^{r) is that for Po('^) both bonds are known to be bonds 
of the same polymer ring while the interaction integral eq [15] corresponds only to a probability 
G{r) for both bonds being in the same chain times a probability G{r) that this chain is closed. 
That Po{r) is negative is of course due to the closure constraint which corresponds to an entropic 
spring force bending the second bond back to the origin. Since this force is scale free (for infinite 
chains) this yields a power law. 
Linear chains. It follows immediately from eq [20] that for linear chains 

P(r) -IM- J_ J_ (22) 
^ ' ~ G{r) ~ 12^ pr^ ^^^^ 

which demonstrates finally the key claim (eq[5]) made in the Introduction (assuming Cqo = 1). 
The normalization factor G{r) is due to the fact that for P(r) both bonds are known to belong 
to the same chain. As compared to the closed cycles the correlation has the opposite sign since 
the attractive spring of the ring closure (indicated by G{qi) in eq [T7|) has been replaced by the 
effective repulsion (indicated by —G{qi)v{qi)G{qi) = —G{qi)/2p in eq[T8]). This repulsion bends 
the second bond away from the origin increasing thus the bond-bond correlation function. 
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Sum rule. Interestingly, the perturbation result, eg 1191 may be rewritten as 

P(r) + ^Po(r) = (23) 

where we have used the normalization factors mentioned above. This "sum rule" suggests a 
geometrical interpretation of the observed relation between infinite linear chains and closed 
cycles which may remain valid beyond the one-loop approximation used here. The idea is that 
in an hypothetical ideal melt containing both linear chains and closed cycles all correlations 
disappear (on distances much smaller than the typical chain sizes) when summed up over the 
contributions of both architectures. The weight {G{r) /2) / p corresponds to the fraction of bond 
pairs in closed loops. Since the orientational correlations in ideal cycles are necessarily long- 
ranged due the ring closure (eg [2T]). it follows, assuming the sum rule, that the same applies to 
bond pairs of linear chains. Since bonds in closed cycles are anti-correlated iPo{r) < 0), they 
must be aligned {P{r) > 0) for linear chains. 

Flory size-distributed chains. We turn now to the upper boundary indicated in eq [5] and 
attempt to characterize finite- effects for u S> 1. We start by considering self-assembled 
chains (branching of chains and formation of closed loops being disallowed) having an annealed 
size-distribution^^ with an exponentially decaying number density = pp'^e~^^ for polymer 
chains of length N with = 1/ {N) being the chemical potential. This so-called "Flory size- 
distribution" is relevant to eguilibrium polymer systems.^' For Flory distributed Gaussian 
chains the form factor becomes"*^^ ^iK) — ^/((fl^)^ + A*) in the intermediate wavevector regime. 
Since the first term in eq[T2]can again be neglected in the incompressible limit (^fc <C 1), this 
yields an effective excluded volume 

i.e. the term (afc)^ in eg [2] is replaced by (afc)^ + p.. Since this applies also for the propagator, 
which becomes G{k) = l/((aA;)2 + /i), the central eqH remains valid for Flory distributed chains 
and, hence, also the sum rule eg 1231 We compute as before eg [20] using now the inverse Fourier 
transform of G{k) in real space 

°<'-' = (^^> 

This yields P{r) = Ii{r)/G{r) = Pa.{r) hi{x) with Pair) denoting the asymptotic power law 
eg [5] (with Cqo = 1) and 

hi{x) = (1 + 2x)%xp(-2x) 

1 /2 

a rapidly decaying function of x = ^yjlr/2a k, r/R with R k, h {N) ' being the typical end-to- 
end distance of the polydisperse system. 

Interestingly, the diagram (1) is not sufficient to characterize the bond-bond correlation for 
larger distances since the last diagram (2) of Figure [U^b) corresponding to the convolution 
integral 

x[-v{k)]G{k). (26) 

provides, as we shall see, the actual cut-off of the power law in this limit. Using again eq [TT] and 
eg [13] the integral factorizes 

Hq) = ^ J J^sGik) . Giq)^ (27) 
= -C2X {aqfG{q) (28) 
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where we have introduced in the last hne the convenient dimensionless constant 

illa)^ f dk o^,,, , , 

« = ii;i3 / (29) 

in which we dump local physics at large wavevector k. Before evaluating the angular correlations 
in real space it is important to clarify the physics described by the diagram. The underlined 
second factor in eq [27] characterizes the alignment of the bond vectors of the monomers hq and 
ni — 1 at a fixed distance r of the monomers no and ni = no + s (Figure [T]J^ a)). Obviously, even 
for Gaussian chains these two bonds become more and more aligned if the distance r = |r| gets 
larger than i.e. when the chain segment becomes stretched. For perfectly Gaussian chains 

the bonds Iq and li at no and ni would still remain uncorrelated, however, since the second 
bond is outside the chain segment on which we have imposed the distance constraint. 

As indicated by the dashed line in the diagram, it is then due to the effective interaction 
between the monomers within the stretched segment (n < ni) and the monomers outside (n > 
ni) that the bonds at ni — 1 and ni get aligned and then in turn the two bonds at no and 
ni. We note that, strictly speaking, C2 depends on the mean chain length (A^), since G{k) is 
a function of /x. However, one checks readily that this effect can be neglected for reasonable 
mean chain lengths. We also note that the constant C2 is finite, since the UV divergence which 
formally arises for large k (where C2 ~ k) may be regularized by local and, hence, model 
dependent physics. We will determine C2 numerically from our simulations of self-assembled 
linear equilibrium polymers (Section (H). 

Assuming a finite and chain length independent coefficient C2 in eq [28] and inserting the prop- 
agator G{q) = l/((ag)^ -|- fi) for Flory distributed chains we obtain by inverse Fourier transfor- 
mation 

hir) = C2 {l^G{r) - 6{r)) (30) 

for the interaction integral in real space. Normalizing l2{'r) as before with G(r) and summing 
over both diagrams this yields 

P(r) =Pa(r)/ii(x)+C2/i (31) 

for r ^ ^ > 0. Comparing both terms in eq [31] one verifies that the crossover occurs at 
r* ~ 6 (N)^^^ in agreement with eq[7] The bond-bond correlation function of an incompressible 
solution of Flory distributed polymers becomes thus constant for r ^ r* . This remarkable result 
is essentially due to the polydispersity. This allows to find for all distances r pairs of bonds Iq 
and li stemming from segments which are slightly stretched by an energy of order /i ^ 1 and 
which are, hence, slightly shorter than a unstretched segment of length s ~ {r/h)"^. Since there 
are more shorter chains and chain segments this just compensates the decay of the weight due to 
the weak stretching. Although the number of such slightly stretched segments decays strongly 
with distance, their relative effect with respect to the typical unstretched segments, e'^ — 1 ~ /i, 
remains constant for all r. It is for this reason that the chemical potential appears in the second 
term of eq 1311 Please note that bond pairs from strongly stretched segments (corresponding to 
an energy much larger than /x) are, however, still exponentially suppressed and can be neglected. 
As we will show now, this is different for monodisperse chains where strongly stretched chain 
segments contribute increasingly to the average for large distances. 

Finite chain size effects: Monodisperse chains. The bond-bond correlation function of 
monodisperse polymer melts may be obtained by inverse Laplace transformation of the result 
obtained for Flory distributed grand-canonical polymers. Note that a very similar calculation 
has been described in detail in ref 13 for the coherent structure factor. The interaction integral 
In{(i) for monodisperse chains of chain length N in reciprocal space may we written 

lN{q) = C-^ [/i-' (/i(g, li) + h{q, /i))] (32) 
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where /(/i) = C[f{N)] = f{N)e~^'^ denotes the Laplace transform of a function f{N). 
The extra factor /i^^ stands for the danghng tails in both diagrams which accounts for the 
combinatorics necessary due to the finite chain length. Ii{q,fi) and 12(9, /x) are the interaction 
integrals computed in the previous paragraph using eqlMlfor the effective interaction potential 
v{k,jj,) and the corresponding Fourier-Laplace transformed propagator G{k,^) = l/((oA;)^ + fi). 
Note that the first term in eq [32] is accurate up to finite-size corrections due to the use of eq [2l] 
for the effective potential. 

We compute then the inverse Fourier transformation /Ar(r) = [Ij\^{q)] and normalize it 
consistently with !F~^ [C~^ [fi~'^G{q, fi)]] using eq[25l This yields 

P{r) = P,{r) hi{x) + ^ h2{x) (33) 

with X = r/{2a^/N) and C2 defined as for Flory distributed chains (eq I29p. The first term due 
to diagram (1) of Figure [Hb) contains a (rapidly decaying) cut-off function 

^ ^ ^ i^erfc(2a;) + 2x i"^erfc(22;) + x^erfc(2x) 

i^erfc(a;) 

with i"erfc(x) denoting the repeated integral of the complementary error function. This func- 
tion is non-monotic and goes through a maximum with an overshoot of about 54% at x ~ 0.39. 
The second term in eq[33]is again due to diagram (2). The function 

/i2(x) = ^^^^^^^ 
4i^erfc(x) 

becomes constant for small x where h2{x) ~ 1 + 2x/-^/7r. We note that as for Flory distributed 
chains the second term in eq 1331 becomes dominant on scales r ^> r*{N) ~ bN^^^. Interestingly, 
for large x we have /i2(x) ~ — 1/2 and P(r) is, hence, predicted to increase as 

P{r) ~ ~ {r/Nf for x > L (34) 

We remind that for a chain segment of arc-length s stretched between its end monomers uq and 
ni one expects (Zq ■ 'i) ~ C2(r/s)^ if r ^ fes^/^. The coefficient C2 stands for the correlation of 
the bond at monomer ni — 1 within the stretched segment with the bond li at the monomer 
Til outside the segment (Figure [T][|b)). The limiting behavior, eq[3U is thus expected since more 
and more bond pairs from strongly stretched chain segments with s ^ N must contribute to 
the average at distances r ^ hN^/'^ . 



Finite persistence length effects. Up to now we have supposed that the chains are perfectly 
flexible. For semifiexible chains (cqo > 1) the above results apply now to the Kuhn segments 
of the chains. The bond length / of the Gaussian reference chain used in the perturbation 
calculation corresponds to the length of the Kuhn segment /k = byjc^o = Icoo (i-e. not to the 
effective bond length 6), the arc- length s to the number of Kuhn segments sk = s/cqo and the 
density p to the density of Kuhn segments 

/O =^ PK = /O/Coo- (35) 

If the bond-bond correlation function P[r) calculated in terms of Kuhn segment units is reex- 
pressed in the natural microscopic units, eq [35] introduces the additional prefactor Cqo indicated 
in eq ([5]). 
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Figure 2: Bond-bond correlation function P(r) obtained for one chain length N = 8192 and 
various monomer densities p as indicated in the panels: (a) Unsealed raw data, (b) rescaled 
correlation function P{r) / Ps,{^) versus the natural scaling variable u = r/^ using the length 
scales ^(/o), h{p) and l{p) indicated in the inset of the panel. The predicted power-law asymptotes 
for dense (eq[5]) and dilute (eg [36j) systems are indicated by bold and dashed lines, respectively. 
The increase of P{r) visible for large distances due to the enhanced weight of stretched chain 
segments is well described by eq[33]as indicated for p = 0.5/8 (thin lines) assuming C2 = 0.14. 

3 Algorithmic and technical issues 

The theoretical predictions derived above are supposed to hold in any dense polymer solution 
containing sufficiently long chains. The numerical data presented below (Figures [2]|31) have 
been obtained using the well-known "bond fluctuation model" (BFM)^"^'^^'^^"^^ — an efficient 
lattice Monte Carlo scheme where a coarse-grained monomer occupies 8 lattice sites on a simple 
cubic lattice (i.e., the volume fraction is 8/9) and bonds between adjacent monomers can vary 
in length and direction. ^"^ All length scales are given in units of the lattice constant and we set 
k-Q = T = 1. Even the partial overlap of monomers is forbidden in the classical formulation of 
the BFM.''' ^^'■^^'^^ We have put the predictions to a test by simulating systems having either a 
quenched and monodisperse or an annealed size distribution: 

(i) Monodisperse systems have been equilibrated and sampled using a mix of local, slithering 
snake, and double bridging Monte Carlo moves. See ref 7 for details. Systems with chain lengths 
up to = 8192 have been obtained for various densities, as indicated in Figure [21 up to the 
"melt density" p = Q.b/SJ We use periodic simulation boxes of linear size L = 512. For the 
largest density used this corresponds to nmon = 2^^ ~ 10'' monomers and to 1024 chains of 
length = 8192 per simulation box. The smallest density indicated {p = 0.00048/8) refers 
to one single chain in the box allowing us to characterize properly the dilute reference point. 
The scaling of bond-bond correlation function with chain length N aX p = 0.5/8 is presented in 
Figure [3l 

(a) As sketched in Figure 01 systems with annealed size distribution — so-called "equilib- 
rium polymers" — have been obtained by attributing a finite and constant scission energy 
E to each bond which has to be paid whenever the bond between two monomers is broken. 
Standard Metropolis Monte Carlo is used to reversibly break and recombine the chains.^' ^^'^^ 
Branching and formation of closed rings are explicitly forbidden. As one expects from stan- 
dard linear aggregation theory, the density of chains pN shows essentially a Flory distribution, 
Pat ~ exp(— A^/ {N)), with the mean chain size (A^) scaling as (A^)^ ~ p exp{E / k^T) }^ Only lo- 
cal hopping moves have been used for sampling equilibrium polymer systems, since the breaking 
and recombination of chains reduce the relaxation times dramatically compared to monodisperse 
systems. 
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4 Simulation results 



Density effects. The scaling of the bond-bond correlation function with density for monodis- 
perse chains is addressed in Figure El Only data for our largest chain length = 8192 is 
presented here to focus first on the discussion of the large-A^ limit. The finite-A^ effects, visible 
nevertheless ("stretching") from the increase of P(r) for large r, will be further discussed below 
at the end of this section (Figures [3] and H]) . As can be seen in panel (a), model-depending 
physics not taken into account by the theory obviously becomes relevant for short distances 
corresponding to segments of a couple of monomers. For clarity we have thus omitted data 
points with r < 5 in all other figures and panels below. 

The power-law behavior observed for small densities is implicit to the swollen-chain statistics 
in dilute good solvents where the root-mean-square segment size R{s) = bos'^^ is known to scale 
with a Flory exponent i^q « 0.588 (with bo denoting the respective power-law amplitude). Since 
P{s) ~ dlR^{s) ~ s~'^" with ujQ = 2 — 2vq (see ref 5) it follows from eqlHthat the bond-bond 
correlation function should scale as 

P{r) « Po(r-) = {bo/lf{r/ho)-^^'^^ ~ r-^'^o (36) 

as indicated by the dashed lines in Figure [21 The index indicates that this is the expected 
asymptotic behavior for asymptotically long chains in the dilute good solvent limit. 

As the density increases, the correlation function still coincides with the dilute behavior for 
distances r smaller than the screening length ^{p) of the semidilute solution where each chain 
segment interacts primarily with itself. (The screening length is indicated in the inset of panel 
(b).) At distances r ^ £,{p)i where the chains overlap and form a "melt of blobs" ,^ the correlation 
function decreases much faster, however not exponentially as one might expect from Flory's 
hypothesis or the (5(r)-scenario ( "self- kicks" ) mentioned in the Introduction, but with a second 
power-law regime with exponent w/i^ = 3. 

That the crossover between both power-law regimes occurs indeed at r ~ ^(p) can better be 
seen from panel (b) where we have replotted the data as a function of the reduced distance 
r/^(p) as suggested by eq[6l Taking apart the finite chain-size effects for large distances we find 
an excellent scaling collapse of the data considering the large spread of the raw data (panel (a)) 
and that no free shift parameter was used. The three independently determined length scales 
^, b and I used for the rescaling according to eqOare compared in the inset of panel (b). The 
latter two lengths are relevant due to the stiffness parameter Cqo = {b{p)/l)'^ which is needed 
(eq[S]) for the vertical scale -Pa(0- Note that the bond length I (triangles) is essentially constant 
for all densities, at least on the logarithmic scales addressed here. The effective bond length b 
approaches / from above, i.e. Coo(p) decreases with increasing p. The dashed line in the inset 
corresponds to the power law expected from the scaling theory of semidilute solutions^ 

b\p) ~ p-{2-o-l)/(3.o-l) ^ p-0.23_ (37) 

Thus, the density dependence of Cqo can not be neglected. The screening length ^ has been 
determined assuming eq [T3] and using the directly measured effective bond length b{p) and 
dimensionless compressibility g{p). We have crosschecked these values with the decay of the 
total structure factor where this has been possible. For not too high densities our data is nicely 
fitted by the power law (bold line) 

^(p) ^ p-o/(3.o-l) ^ p-0.77 (38) 

as expected for semidilute solutions.^ Obviously, the semidilute power-law relations cannot hold 
strictly for the highest densities where the compressibility (and, hence, the blobs size) becomes 
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Figure 3: P{r) for r > 5 for different chain lengths at melt density p = 0.5/8. The thin lines 
indicate eq[33]for each N assuming C2 = 0.14. The limiting behavior for large distances, eq[3ll is 
given by the dashed lines, (a) The data approach systematically the power-law decay (bold line) 
predicted by eqOwith increasing chain length A^. (b) P{r)N as a function of r/r* assuming 
eqO (main panel) and as a function of x = r /{laN^^"^) (inset). The first scaling variable scales 
the data around the minimum of P{r), the second for larger distances. 



too small. "'^ However, the differences are small on logarithmic scales and we obtain a very 
similar scaling plot by insisting on eqs [37] and [38] for all densities (not shown) . 

The density crossover scaling implies obviously the matching of the dilute and dense asymptotic 
power-law predictions 

PoiO^m)- (39) 

As can be checked using the well-known scaling relations 6^ ~ C^/fi'; C'^P ^ 9 and ^ ~ bog'^" 
defining the semidilute solution^ (which are also implicit to eq [37] and eq [38]) , eq [39] requires 
the prefactor Coo{p) in eq U We have checked that Po(0 ~ scales the data, while Pa (6 

without Coo does not (not shown). The successful scaling collapse confirms, hence, the rescaling 
of the Kuhn segments presented at the end of section [5] 

Finite A^-efFects for monodisperse chains. We concentrate in the reminder of this section 
on dense melts with p = 0.5/8. Chain length effects are investigated in Figure [3] for monodis- 
perse chains. The unsealed raw data are presented in panel (a). As expected from theory, P{r) 
increases for large distances due to the increasing weight of stretched segments contributing to 
the average at distance r. With increasing chain length A^ this effect becomes less important, 
however, and our data approaches systematically the asymptotic power-law decay (eq[5|) indi- 
cated by the bold line. The thin lines represent the complete prediction, eg 1331 for different A'^. 
Numerical data and theory agree nicely, especially for large chains. The deviations observed 
for A^ < 256 and for fully stretched chains (r ~ IN) are, of course, expected. Please note that 
the cut-off function hi{x) of the asymptotic power law obtained by diagram (1) may be ignored 
{hi{x) — > 1) without changing the plot. The finite- A^ effects are in fact dominated by the func- 
tion h2{x) of the second term in eq[331 We remind that eq [M] predicts ultimately a quadratic 
increase with distance of the bond-bond correlation function as indicated by the dashed line for 
one chain length (A'^ = 64). It is essentially the limited simulation box size (L = 512) of the 
present study preventing us from demonstrating numerically this asymptotic behavior which 
should be accessible otherwise for the longer chains (A'' S> 10^). 

The scaling with chain length is addressed in panel (b). Obviously, eq[33]is not compatible with 
one scaling variable allowing to collapse all data. The crossover from the power-law asymptote 
(diagram (1)) to the A^-dependent correction (due to diagram (2)) is, however, well described 
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Figure 4: Bond-bond correlation function P{r) for equilibrium polymers at density p = 0.5/8 of 
the BFM for different bond energies E. This energy has to be paid for the (reversible) breaking 
of a bond as indicated by the dashed line in the sketch where a chain of = 9 monomers is 
broken into two chains of length = 4 and N = 5. The corresponding mean chain lengths 
{N(E)) are indicated. If P{r) (N) is traced as a function of r/r* with r* = b{N)^^^ all data 
points collapse. Note that P{r) {N) — > C2 = 0.14 for r/r* ^ 1. The complete prediction eg [311 
(thin line) interpolates perfectly between the power-law asymptote (bold line) and the plateau 
(dashed line). 



by plotting P{r)N as a function of r/r* with r* = bN^^^ using as a scale the minimum of eg 1331 
As shown in the main panel this yields a numerically satisfactory scaling over two orders in 
magnitude of the reduced distance r/r*, especially for our largest chains. This scaling fails for 
large distances (r ^ R{N)), however, as emphasized by the theoretical predictions (thin lines) 
indicated for A^ = 64 and A^ = 8192 (bottom). In this limit only the second term in eg [331 
matters and, as can be seen from the inset in Figure [3l^b) , a data collapse can be achieved by 
tracing P{r)N as a function of x ~ r/R{N). As already noted, the observed deviations from the 
limiting behavior for x ^ 1 (eg I34p are expected (i) due to the breakdown of the Gaussian chain 
model at r ~ A^^ for small chain lengths and (ii) due to the restricted box size of the present 
study. 

Equilibrium polymers. Fortunately, this scaling issue is much simpler for equilibrium poly- 
mers as shown in Figure [3] where we have plotted P{r) (N) as a function of r/r* with r* = 
b (N)^^^ using the indicated mean chain lengths (A^). Note that the error bars (not shown) 
become clearly much larger than the symbol size for large bond energies E > 12. It is fair to 
state, however, that all data points collapse nicely on the one master curve indicated by the 
thin line obtained from eg [31] for (A^) — > oo. That the variable x effectively drops out stems 
from (i) the rapid decay of the exponential cut-off function of the first term in eg [31] and (ii) the 
^-independence of the correction term. That the bond-bond correlation function of equilibrium 
polymers becomes constant for large distances confirms a non-trivial prediction of the theory. 
We used the clearly visible plateau to determine the coefficient C2 = 0.14 required for the in- 
teraction diagram (2) and already used in the previous Figures [2] and [31 We note finally that 
this best fit value of C2 is rather close to P{s = 1) = {In ■ In+i) /^^ ~ 0.10, the independently 
determined bond-bond correlation between adjacent bond vectors. 
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5 Conclusion 



Summary. Focusing on dense solutions of linear and (essentially) flexible polymer chains we 
have investigated analytically and by means of Monte Carlo simulation the intrachain angular 
correlations with respect to the distance r between bond pairs. Motivated by recent work showing 
the power-lay decay of the bond-bond correlation function P{s) with curvilinear distance s (eq[TJ 
we addressed the question whether the correlations are indeed long-ranged in space or only due 
to the return probability ("self-kicks") of a chain. Our calculations of P[r) have been based 
on a standard one- loop perturbation scheme (Figure Hb)). The power-law decay predicted 
for asymptotically long chains (eq [5|) is well confirmed by our numerically data approaching 
systematically the asymptotic envelope Ps.{r) with increasing density (Figure [2|) and (mean) 
chain length (Figures [3] and UD . As postulated in eq El density effects are found to scale as 
P(r)/Pa(0 = f{u = r/^) with ^{p) being the independently measured screening length. For 
u <C 1 we confirm the expected power law for dilute good solvents (eq l36p . More importantly, the 
explicit compressibility dependence of the bond-bond correlation function drops out for u ^ 1 
where f{u) =^ i.e. polymer solutions behave on large scales as incompressible packings of 

blobs and this for all densities provides that the chains are sufficiently long. Finite-chain size 
effects are also successfully predicted by the theory for both monodisperse polymers (eq [33]) and 
equilibrium polymers (eq [3T]). It should be emphasized that the fit to the asymptotic power law 
is parameter free and that fitting the A^-effect only requires one additional free parameter, C2 
(eq [28]). due to local orientational correlations between adjacent bonds which regularize diagram 
(2) in Figure mb). 

Outlook. Interestingly, the presented perturbation calculation for dense polymer chains may 
also be of relevance to the angular correlations of dilute polymer chains at and around the 
0-point which has received attention recently. The reason for this connection is that (taken 
apart different prefactors) the same effective interaction potential (eq[5]) enters the perturbation 
calculation in the low wavevector limit. We expect therefore similar genuinely long-ranged 
angular correlations for asymptotically long 0-chains (as implied by eq[3D rather than the "self- 
kicks" suggested in the literature. Strong finite-A effects are again expected, however, and 
much longer chains as the ones presented in the current numerical studies of 0-chains are required 
to demonstrate the asymptotic power-law behavior. 

The present study has focused on the first Legendre polynomial of the intrachain bond-bond 
correlations. Since this is not an easily experimentally accessible property it should be men- 
tioned that higher Legendre polynomials have also been predicted following similar perturbation 
calculations as the ones presented above. For instance it is possible to show that the second 
Legendre polynomial should decay as the fifth power of distance if averaged over all intra-chain 
contributions and as the sixth power if averaged over all bond pairs at a given distance. Un- 
fortunately, we are at present unable to demonstrate these predictions numerically due to the 
stronger power-law decay requiring much better statistics. This study is currently underway. 

Our one-loop perturbation calculation show that for infinite chains or Flory distributed (grand- 
canonical) polymers the bond-bond correlations observed in systems of linear chains are equiv- 
alent to the subtraction of angular correlations due to closed cycles (eq [23]) . The same sum 
rule can be shown to hold for higher Legendre polynomials summing over intra- and inter-chain 
contributions. This finding suggests a deeper purely geometrical description of the observed 
long-ranged orientational correlations relating this issue to the recently discovered Anti-Casimir 
forces in polymer melts^^'^^ arising due to a similar subtraction of soft fiuctuation modes, not 
present in the linear polymer system but in its hypothetical counterpart containing both chain 
architectures. Taking advantage of the polymer-magnetic analogy^ we plan to address this 
fundamental connection in a more theoretical paper. 
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